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ABSTRACT 

The Herschel Multi-tiered Extragalactic Survey, HerMES, is a legacy program designed 
to map a set of nested fields totalling ~ 380 deg 2 . Fields range in size from 0.01 to 
~ 20 deg 2 , using Herschel-SPIKE (at 250, 350 and 500 /xm), and Herschel - PACS (at 
100 and 160 /mi), with an additional wider component of 270 deg 2 with SPIRE alone. 
These bands cover the peak of the redshifted thermal spectral energy distribution from 
interstellar dust and thus capture the re-processed optical and ultra-violet radiation 
from star formation that has been absorbed by dust, and are critical for forming a 
complete multi-wavelength understanding of galaxy formation and evolution. 

The survey will detect of order 100,000 galaxies at 5cr in some of the best studied 
fields in the sky. Additionally, HerMES is closely coordinated with the PACS Evolu- 
tionary Probe survey. Making maximum use of the full spectrum of ancillary data, 
from radio to X-ray wavelengths, it is designed to: facilitate redshift determination; 
rapidly identify unusual objects; and understand the relationships between thermal 
emission from dust and other processes. Scientific questions HerMES will be used to 
answer include: the total infrared emission of galaxies; the evolution of the luminosity 
function; the clustering properties of dusty galaxies; and the properties of popula- 
tions of galaxies which lie below the confusion limit through lensing and statistical 
techniques. 

This paper defines the survey observations and data products, outlines the primary 
scientific goals of the HerMES team, and reviews some of the early results. 

Key words: surveys - infrared: galaxies - submillimetre: galaxies - galaxies: evolu- 
tion 


1 INTRODUCTION & SCIENCE GOALS 

Understanding how galaxies form and evolve over cosmolog- 
ical time is a key goal in astrophysics. Over the last decade 
our understanding of the background cosmology has im- 
proved to such an extent (e.g. Spergel et al. 2003) that we 
think we have a reasonable understanding of the formation 
of structure in the underlying dark matter distribution (e.g. 
Springel et al. 2006). However, galaxy formation and evo- 
lution are driven by dissipative, non-linear processes within 
the potential wells of virialized dark matter halos which are 
much more complex physical processes which have defied full 
modelling. Observations play a critical role in constraining 
models of galaxy formation, the evolution of star-formation 
activity, and the various roles played by galaxy stellar mass, 
dark matter halo mass, and environment. 

The central importance of far-infrared (FIR) and sub- 
millimetre (sub-mm) observations becomes clear when one 
realizes that the approximately half of all the luminous 
power (Puget et al. 1996; Fixsen et al. 1998; Lagache et al. 
1999) which makes up the extra-galactic background radia- 
tion - power which originated from stars and active galac- 
tic nuclei (AGN) - was emitted at optical/ultraviolet wave- 
lengths, absorbed by dust, and reradiated in the FIR/sub- 
mm. To form a complete picture of the evolution of galaxies, 
the optical regime alone cannot be used to fully trace the 
activity (e.g., the brightest sub-mm galaxy in the Hubble 
Deep Field is not even detected in the optical Dunlop et al. 
2004). Furthermore, sub-mm observations provide a bridge 
in both wavelength and redshift between the z > 2 Universe, 
primarily probed on the Rayleigh- Jeans side of the spectral 
energy distribution (SED) by ground based sub-mm tele- 
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Figure 1. Model spiral (green), star-burst galaxy (blue) and 
AGN (red) Spectral Energy Distributions (SEDs) normalised to 
the same Lfir and plotted in their rest-frame with the Herschel- 
PACS and Rersc/ieZ-SPIRE bands at A = 100, 160, 250, 350 and 
500 /im plotted at A/(l + z) for a galaxy at z = 1.5. Note that 
the He rs ch e l- S P THE band at 250 //remeasures a similar flux den- 
sity for all and so is a reasonable proxy for the Lfir for these 
templates. 


scopes, and the lower- z Universe, sampled on the Wein side 
of the SED by Spitzer. 

FIR/sub-mm luminosity is thought to arise primarily 
from dust heated by the massive stars in star formation re- 
gions and so may be used as a direct estimate of star forma- 
tion activity. Additional contributions are expected to arise 
from dusty tori surrounding AGN at shorter wavelengths, 
and there may be non-negligible contributions from the il- 
lumination of dust by evolved stars. 

Previous surveys from space-based observatories: IRAS 
(e.g. Saunders 1990; Oliver et al. 1992); ISO (e.g. Elbaz 
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et al. 1999; Oliver et al. 2002, and references therein); and 
Spitzer (e.g. Shupe et al. 2008; Frayer et al. 2009, and refer- 
ences therein); and at sub- mm wavelengths from the ground 
with: SCUBA at 850 /im (e.g. Eales et al. 1999; Hughes et al. 
1998; Smail et al. 1997; Coppin et al. 2006), Bolocam (Mal- 
oney et al. 2005a, e.g.); SHARCII (e.g. Khan et al. 2007); 
MAMBO (e.g. Greve et al. 2008); LABOCA (e.g. Weifi et al. 

2009) ; and AzTEC (e.g. Scott et al. 2010b), demonstrated 
strong evolution in galaxies at both mid-infrared (MIR) and 
FIR wavelengths. This evolution is attributed to a decline 
in the average star-formation density with time, and partic- 
ularly a decline in the role of the more luminous infrared 
galaxies (LIRGs), which are thought to be the progenitors 
of massive galaxies today (e.g. Le Floc’h et al. 2005). 

This strong evolution has been challenging for physical 
models of galaxy formation to reproduce. They find they 
must invoke drastic modifications, such as altering the ini- 
tial mass function (e.g. Baugh et al. 2005), in order to match 
these observations as well as optical and near infrared con- 
straints on the stellar mass today. 

Using a different approach, phenomenological galaxy 
population models attempt to describe what is currently ob- 
served and also predict what we would expect for Herschel. 
Different groups use different combinations of galaxy popu- 
lations to reproduce the observations; for example, Lagache 
et al. (2003, and Fig. 2) use two peaks of luminosity density 
at z ~ 1 and z ~ 2 to describe the data, which are not seen 
in other models. Such differences between the pr e-Herschel 
models indicate the lack of constraint on the spectral energy 
distributions and redshift distributions. 

The potential of sub-mm surveys has been demon- 
strated by the BLAST telescope (Devlin et al. 2009). BLAST 
was a balloon-borne telescope with a focal plane instru- 
ment based on the SPIRE (Griffin et al. 2010) photometer 
design and using similar detectors tailored to higher pho- 
ton loading, and was a successful technical and scientific 
pathfinder for SPIRE on Herschel, probing the wavelength 
regime where the SED of redshifted galaxies and the infrared 
background peak. 

The Herschel Space Observatory (Pilbratt et al. 2010) is 
carrying out surveys of unprecedented size and depth, vastly 
improving the state of observations in this under-explored 
waveband. The imaging instruments SPIRE (Griffin et al. 

2010) and PACS (Poglitsch et al. 2010), which together fully 
constrain the peak of the FIR/sub-mm background, allow us 
to thoroughly investigate the sources in the infrared back- 
ground and characterize their total obscured emission (see 

e-g- Fig- !)• 

The Herschel Multi-tiered Extra-galactic Survey (Her- 
MES 1 ) is a Guareenteed Time Key Program on Herschel 
which will provide a legacy survey of star forming galax- 
ies over the wavelengths at which the galaxies and infrared 
background peak. The majority of science goals require 
multi-wavelength support and the fields we have chosen are 
among the best in the sky for multi-wavelength coverage (see 
Section 4.2) maximising their legacy value. 

In Section 2 we define the survey. In Section 3 we de- 
scribed some of our goals and early results. In Section 4 we 


1 http://hermes.sussex.ac.uk. Hermes is also the Olympian mes- 
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outline our expected data products and delivery time-scales 
before concluding in Section 5. 


2 SURVEY DESIGN 

Our survey is defined by Astronomical Observing Requests 
(AORs). For convenience we have grouped the AORs by 
sets, which in this paper are identified with numbers, e.g., 
ELAIS N2 SWIRE is #41. A summary of the AOR sets 
is given in Table 1. Details of the observing modes can be 
found in the Herschel observers’ manuals (available from 
http : //herschel . esac . esa. int /Document at ion. shtml). 

Detector hits maps 2 , which accurately define the cover- 
age of the survey and should be used for any detailed plan- 
ning of complementary surveys, are provided on our web site 
http://hermes.sussex.ac.uk. We also provide files which 
define the approximate boundaries of homogenous regions 
(e.g. as marked in Fig. 5). These survey definition products 
are updated as the survey progresses. Our sensitivities have 
been quoted using official mission values given in Table 3. 

The current AORs which define our program 
can be retrieved from the Herschel Science Archive 
http : //herschel . esac . esa. int/Science_Ar chive . shtml 
using HSpot and the proposal IDs SDP_soliver_3 and 
I<PGT_soliver_l and GT2_mviero_l. 

Here we summarise the basis of our survey design. 

2.1 Requirements 

HerMES was designed to fulfil multiple objectives, which 
are outlined in Section 3. The Herschel bands can probe the 
peak of the far infrared spectral energy distributions of star 
forming galaxies and thus measure the infrared luminosity, 
Lir, see Figure 1 and Table 2. Our primary criterion was 
to sample the (Lir, z) plane of star-forming galaxies uni- 
formly and with sufficient statistics to a redshift of 0 < z < 3. 
Specifically, we take a bin resolution of AlogLiR Az = 0.1 
(e.g. AlogLiR = 0.5, Az = 0.2) and require 75 galaxies per 
bin to give 12 per cent accuracy (or 10 per cent accuracy 
when further divided into three sub-samples). This resolu- 
tion corresponds to the scale of features in the luminosity 
density surface from the Lagache et al. (2003) model, for 
example. Using the model luminosity functions we can cal- 
culate the area needed to reach this goal for each luminos- 
ity and redshift. Each tier thus probes a given (Lir, z) re- 
gion bounded by the areal constraint and the flux limit (see 
Fig. 2). An optimized sampling over wavelength is achieved 
by combining HerMES with the PACS Evolutionary Probe 
survey (PEP, Lutz et al. 2011). 

HerMES was thus designed to comprise a number of 
tiers of different depths and areas (Tables 5 and 7). Her- 
MES samples the higher luminosity objects, which are bright 

2 These maps and Table 1 gives coverage for SPIRE observations 
as counts of 250 pm detector samples per 6”x6” pixel. This can 
be converted to a bolometer “exposure” time per pixel by dividing 
by the sampling frequencies of 18.6 Hz for SPIRE scanning at 
nominal and fast rates and 10 Hz for parallel mode. The hits 
in other arrays can be estimated by scaling by the numbers of 
detectors in the arrays (129, 88, 43) and the pixel sizes (6", 10”, 
12 "). 
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Figure 2. Far infrared luminosity density in log 10 (LQ/i — 3 Mpc 3 dex — - 1 ) (grey-scale and contour diagram) as a function of far infrared 
luminosity (rr-axis) and redshift (y-axies) - from the model of Lagache et al. 2003. The power of different survey elements to probe this 
space are indicated by overlays. Each panel shows survey elements at different wavelengths; reading left-to-right from the top they are 
100, 160, 250, 350 and 500 fim. Surveys are deemed to properly sample the space if they can detect galaxies of these FIR luminosities 
at the 5-cr instrumental noise level and with more than 75 galaxies in bins of AlogL Az = 0.1. These two constraints are marked with 
dotted lines and are hatched. The different survey levels defined in Table 7 are shown with: Levels 2-4 - blue; Level 5 - red; Level 6 
- magenta and HeLMS- - green. Level- 1 (cyan) does not have enough volume to satisfy the number of galaxies criterion and so only 
the instrumental noise limit is shown. The 5 cr confusion noise levels (after 5 a clipping) from Berta et al. (2011, at 100 and 160^m) and 
Nguyen et al. (2010, at 250, 350 and 500/zm) with yellow/black lines. Note the bimodal peaks at z ~ 1 and £ ~ 2.5 
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Set 

Level 

Target 

Mode 

JVaor 

T 

[hr] 

A rep 

Asamp 

h 

['] 

h 

n 

e 

[deg] 

^nom 

[deg 2 ] 

fimax 

[deg 2 ] 

^good 

[deg 2 ] 

DR 

1 

CD 

Abell 2218 

Sp. Nom. 

2 

9.29 

100 

1118 

4 

4 

84 


0.14 

0.10 

SDP 

2 

CD 

Abell 1689 

Sp. Nom. 

8 

1.97 

48 

235 

4 

4 

18 


0.11 

0.08 


3 

CD 

MS0451. 6-0305 

Sp. Nom. 

8 

1.97 

48 

235 

4 

4 

5 


0.11 

0.08 

DR1 

4 

CS 

RXJ13475-1145 

Sp. Nom. 

8 

1.97 

48 

234 

4 

4 

17 


0.11 

0.08 


5 

CS 

Abell 1835 

Sp. Nom. 

8 

1.97 

48 

236 

4 

4 

16 


0.11 

0.08 


6 

CS 

Abell 2390 

Sp. Nom. 

8 

1.97 

48 

235 

4 

4 

81 


0.11 

0.08 


7 

CS 

Abell 2219 

Sp. Nom. 

8 

1.97 

48 

234 

4 

4 

66 


0.11 

0.08 

DR1 

8 

CS 

Abell 370 

Sp. Nom. 

8 

1.97 

48 

233 

4 

4 

70 


0.11 

0.08 


9 

CS 

MS1358+62 

Sp. Nom. 

8 

1.97 

48 

235 

4 

4 

76 


0.11 

0.08 


10 

CS 

C10024+16 

Sp. Nom. 

8 

1.97 

48 

235 

4 

4 

61 


0.11 

0.08 


11 

CH 

MS1054. 4-0321 

Sp. Nom. 

8 

2.18 

16 

131 

15 

10 

22 


0.24 

0.16 


12 

CH 

RXJ0152. 7-1357 

Sp. Nom. 

8 

2.18 

16 

131 

15 

10 

165 


0.24 

0.16 


13 

LI 

GOODS-S 

Sp. Nom. 

76 

20.22 

76 

730 

20 

20 

14 


0.51 

0.35 


22 

L2 

COSMOS 

Sp. Nom. 

24 

50.13 

24 

388 

85 

85 

70 


3.49 

2.82 


14 

L2 

GOODS-N 

Sp. Nom. 

1 

13.51 

30 

416 

30 

30 

42 


0.64 

0.55 

SDP 

15 

L2 

ECDFS 

Sp. Nom. 

19 

8.78 

19 

232 

30 

30 

44 


0.79 

0.58 

DR1 

17 

L3 

Groth Strip 

Sp. Nom. 

7 

3.54 

7 

85 

67 

10 

130 


0.82 

0.60 

DR1 

18 

L3 

Lockman-East ROSAT 

Sp. Nom. 

7 

3.2 

7 

87 

30 

30 

77 


0.77 

0.57 


18B 

L3 

Lockman-East Spitzer 

Sp. Nom. 

4 

4.53 

4 

32 

80 

40 

149 


1.78 

1.40 


19 

L3 

Lockman-North 

Sp. Nom. 

1 

3.91 

7 

104 

35 

35 

1 


0.74 

0.65 

SDP 

23 

L4 

UDS 

Sp. Nom. 

7 

10.54 

7 

110 

66 

66 

20 


2.46 

2.02 


24 

L4 

VVDS 

Sp. Nom. 

7 

10.39 

7 

110 

66 

66 

21 


2.46 

2.02 


22B 

L5 

COSMOS HerMES 

Sp. Nom. 

8 

25.20 

8 

128 

110 

110 

70 


5.04 

4.38 


27 

L5 

CDFS SWIRE 

Sp. Fast 

10 

41.72 

20 

81 

190 

150 

99 


12.18 

11.39 


28 

L5 

Lockman SWIRE 

Sp. Fast 

2 

13.51 

2 

16 

218 

218 

2 


18.2 

17.37 

SDP 

28B 

L5 

Lockman SWIRE 

Sp. Fast 

8 

41.26 

8 

58 

220 

180 

50 


15.26 

7.63 


42 

L7 

HeLMS 

Sp. Fast 

11(10) 

103.4 

2 


1560 

750 

15 

270 




20 

L3 

Lockman-North 

PACS 

12 

13.96 

11 


30 

30 

42 

0.25 



SDP 

20B 

L3 

Lockman-North 

PACS 

20 

20.89 

20 


30 

30 

42 

0.25 




21 

L3 

UDS HerMES 

PACS 

25 

25.93 

14 


30 

30 

0 

0.25 




25 

L4 

UDS 

PACS 

12 

40.19 

7 


57 

57 

0 

0.9 




29 

L5 

EGS HerMES 

Parallel 

7 

22.68 

7 

93 

150 

40 

131 


3.50 

2.67 


30 

L5 

Bootes HerMES 

Parallel 

5 

20.33 

5 

70 

80 

80 

0 


4.21 

3.25 

DR1 

31 

L5 

ELAIS N1 HerMES 

Parallel 

5 

20.82 

5 

72 

95 

95 

38 


3.74 

3.25 

DR1 

32 

L5 

XMM VIDEOl 

Parallel 

4 

13.44 

4 

65 

106 

75 

107 


3.20 

2.72 


32B 

L5 

XMM VIDE02 

Parallel 

4 

8.88 

4 

53 

106 

44 

107 


2.12 

1.74 


32C 

L5 

XMM VIDE03 

Parallel 

4 

13.44 

4 

53 

106 

75 

107 


3.19 

2.73 


33 

L5 

CDFS SWIRE 

Parallel 

4 

50.42 

4 

57 

204 

170 

101 


11.87 

10.89 


34 

L5 

Lockman SWIRE 

Parallel 

4(2) 

71.22 

4 


215 

215 

154 


17.86 

16.08 


39B 

L5 

ELAIS SI VIDEO 

Parallel 

4 

17.72 

4 

56 

138 

80 

87 


4.42 

3.72 


35 

L6 

ELAIS N1 SWIRE 

Parallel 

2 

28.0 

2 

28 

207 

192 

55 


13.37 

12.28 


36 

L6 

XMM-LSS SWIRE 

Parallel 

6 

45.58 

2 

29 

180 

180 

82 


21.62 

18.87 

DR1 

37 

L6 

Bootes NDWFS 

Parallel 

4 

27.99 

2 

30 

243 

80 

145 


11.3 

10.57 

DR1 

38 

L6 

ADFS 

Parallel 

2 

18.11 

2 

28 

190 

122 

80 


8.58 

7.47 

DR.1 

39 

L6 

ELAIS SI SWIRE 

Parallel 

2 

17.9 

2 

28 

140 

81 

91 


8.63 

7.86 


40 

L6 

FLS 

Parallel 

2 

17.1 

2 

29 

160 

138 

5 


7.31 

6.71 

SDP 

41 

L6 

ELAIS N2 SWIRE 

Parallel 

2 

17.1 

2 

26 

177 

119 

147 


9.06 

7.80 



Table 1. Summary of the HerMES observations. The full set of Astronomical Observation Requests (AORs) are available through ESA’s 
Herschel Archive. We have grouped Aaor observations of the same field at the same level made with the same mode and areal size into 
a ‘set’ (the number of AORs still to be scheduled after 2011 Dec 21 is indicated in parentheses). The first five columns in the Table give: 
the set identification number; the design level; the target name, the Herschel observing mode and the number of AORs in the set. T is 
the time used or allocated for this set. A rep is the total number of repeats of the observing mode in the set. All our SPIRE nominal 
(30 // s -1 ) and fast mode (60 // s — 1 ) (Sp. Nom. and Sp. Fast) observations include a scan in the nominal and orthogonal direction, so 
1 repeat is 2 scans). For SPIRE observations that have been executed A samp is the median number of bolometer samples per pixel in 
the 250 fi m map (6 X/ x 6 ;/ pixels). This can be converted to exposure time per pixel or to other bands as described in footnote 2. The 
error per pixel in our SPIRE maps as processed by the standard HIPE pipeline are cr| 50 = cr^/Ar samp with ctq = 896 d= 11, 1554 =b 27 
and ~ 1440 mjy 2 beam - 2 for Parallel, Sp. Nom. and Sp. Fast modes respectively. Ii,l 2 are sides of a rectangle with near homogenous 
coverage. 6 is the roll angle with short-axis of that rectangle measured East of North. For SPIRE observations that have been executed 
fimax is the total area of pixels with any 250 pm coverage and fl goo d is the area of pixels where the number of bolometer samples per 
pixel in the 250 pm map is greater than A sa mp/2. For PACS fields or unobserved SPIRE fields fl n om gives the nominal area of region. 
The final column indicates which observations are included in our data releases; observations marked SDP were released in our Early 
Data Release, observations marked SDP or DR1 will be released in our First Data Release. Set numbers #16 and #26 were removed 
from the programme. 
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at Nominal Wavelength [/rm] 



100 

160 

250 

350 

500 

Instrument 

PACS 

PACS 

SPIRE 

SPIRE 

SPIRE 

Filter name 

Blue2 

Red 

PSW 

PMW 

PLW 

Min A [/xm] 

85 

125 

210 

300 

410 

Max A [/xm] 

125 

210 

290 

400 

610 


Table 2. Basic band information for the different Herschel chan- 
nels used by HerMES. Data is taken from SPIRE and PACS Ob- 
servers’ Manuals V2.4/V2.3 (respectively). 


but rare, in the wide shallow tiers, and the lower luminosity 
galaxies, which are faint but common and confused, in the 
deep narrow tiers. Our design has evolved during the mission 
but since our initial design had cluster observations (nomi- 
nally deep, shallow and high- 2 ) and six nominal levels from 
deep and narrow Level 1 to wide and shallow Level 6 and 
we will maintain those descriptions even though the depths 
have changed. 

Confusion is a serious issue for Herschel and SPIRE 
in particular, and is an important driver in deciding sur- 
vey depth (Table 5). To estimate the confusion level we as- 
sembled galaxy models (e.g. Lagache et al. 2003), compared 
them to existing survey data, and calculated the confusion 
limit using the criteria for source density of 30 beams per 
source and width of the sky intensity distribution. We em- 
ploy a number of techniques to overcome the problem of 
confusion. It is those analyses which motivate the deepest 
tiers: the lensed clusters fields; and the fast scanned elements 
of the wide Level 5 tier. 

An additional consideration is the volume of the sur- 
vey needed for a representative sample of the Universe, to 
provide a sufficient range of environments, and enough inde- 
pendent regions to study clustering (e.g. Fig. 12). Examina- 
tion of each of those requirements requires survey co-moving 
volumes of 10 6 — 10'Mpc 3 or larger. E.g. the number den- 
sity today of halos of dark matter mass Mdm > 10 ls Mg 
is around 10 -6 /i 3 Mpc _3 (Mo & White 2002). This is iden- 
tical to the co-moving number density of their progenitors 
i.e. ~ 0.3 — 0.4 deg -2 for survey shells of A z = 0.1. This 
provides additional motivation for fields of order 10s deg 2 
to provide statistical samples. Sampling variance would still 
be an issue if the smaller deeper levels were contiguous so 
we split these into a number of fields to enable us to reduce 
and estimate the sampling variance errors. 

The SPIRE and PACS depths for the cluster observa- 
tions were determined by the desire to ensure the detection 
of z ~ 3 sources in both the SPIRE 250 and PACS 100 /im 
channels. 

2.2 Choice of Fields 

In order to pursue multi-wavelength analyses, we have se- 
lected fields (Fig. 3, Table 5) which are among the most 
intensively observed at all wavelengths. These incude: radio 
(VLA, WRST, GMRT, ATCA); sub-mm (SCUBA, Bolo- 
cam, AzTEC, MAMBO); mid and far infrared ( Spitzer , ISO, 
AKARI)- near-infrared (UKIRT, VISTA); optical ( HST , 
Subaru SuprimeCAM, CFHT MegaCAM, KPNO Mosaicl, 
CTIO MOSAIC2, INT WFC); UV ( GALEX) and X-ray 
(XMM-Newton, Chandra) . A description of ancillary data is 
given in Section 4.2). Extensive redshift and/or photometric 


redshift surveys are either available or underway for most of 
these fields. 

An additional consideration was that the contamination 
from Galactic emission (or cirrus) should be minimal. The 
larger mirror means that this cirrus is less of a concern for 
extra-galactic surveys with Herschel than it was for Spitzer, 
as discussed in Oliver (2001). This means that our require- 
ment for low-levels of cirrus are automatically satisfied by 
our criteria of good ancillary data, as illustrated in Fig. 3. 

The defining criterion was coverage at mid/far infrared 
wavelengths not accessible to Herschel, or where Herschel 
is relatively inefficient due to its warm mirror. Specifically 
we required Spitzer MIPS coverage at 24 and 70 pm. At the 
time of design the one exception to this was the AKARI 
Deep Field South, which did not have Spitzer coverage but 
did have coverage at 65, 90, 140, and 160 pm from AKARI 
(Matsuura et al. 2010). However, this field has since been 
observed by Spitzer MIPS (Scott et al. 2010a). The HcLMS 
field, which was added in 2011 for studying large-scale struc- 
ture and the bright end of the number counts, does not have 
ancillary data from Spitzer. However, being located on the 
SDSS Stripe 82 region, HeLMS does have ancillary coverage 
from many other facilities. 

A detailed discussion of the specific observations which 
were considered in the design of the fields is given in Ap- 
pendix A 

The deep and shallow cluster targets are well-studied 
strong lenses at modest redshift. They were selected in con- 
sultation with the PEP team - with HerMES carrying out 
the SPIRE observations and PEP the corresponding PACS 
observations. The high-z clusters were selected for environ- 
mental studies also in consultation with the PEP team. 


2.3 Observing Modes 

The mapping of Levels 1-4 (#1,11-19, 22, 23) is performed 
using SPIRE ‘Large Map’ mode. This mode is described in 
detail in the SPIRE Observers’ Manual. 3 This is the default 
SPIRE observing mode for any field size larger than 4'x4'. 
In this mode maps are made by scanning the telescope be- 
cause it eliminates off-beam confusion, aliows measurement 
of extended emission, and increases observing efficiency for 
larger fields. Since our smallest blank field to be mapped 
(Level 1) is ~16 / xl6 ; , this mode was the natural choice for 
our program. 

The SPIRE cluster observations were originally de- 
signed using the ‘Large Map’ mode covering a nominal field 
of 4' x 4' as this was the smallest map that could be made us- 
ing scanning. Abell 2218, #1, was carried out in that mode. 
We moved to ‘Small Map’ mode (#2-10) in which the map is 
made by two short cross-scans with the telescope once that 
became available, as that was more efficient for smail fieids. 

When building maps the telescope is scanned at an an- 
gle of 42.4° with respect to the Z axis of the arrays, (see Fig- 
ure 3.1 and 3.3 of the SPIRE Observers’ manual, V2.4). This 
produces a fully-sampled map, despite the focal plane not 
being fully sampled. The offset between successive scans (or 

3 The SPIRE Observers Manual is avail- 
able from the Herschel Science Centre 

http:/ /herschel. esac.esa.int/Docs/SPIRE/html/spire_om.html 
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Figure 3. Map of dust emission from the Galaxy, with HerMES fields over-plotted. The image is the 100 pm, COBE- normalised, IRAS 
map of extended emission (Schlegel et al. 1998). The projection is Hammer-Aitoff in Galactic coordinates. The sky brightness is plotted 
on a false-colour logarithmic scale, with regions of very low Galactic emission appearing black and the Galactic plane yellow. In addition 
to the blank fields marked, HerMES has also observed 12 known clusters. 



Figure 4. Maps of the number of bolometer samples per pixel of four deep SPIRE 250 pm observations. From left to right: Abell 2218 
which was observed in SDP without dithering; Abell 2219 which was taken with dithering; GOODS-N (taken in SDP without dithering) 
and ECDFS with dithering. FITS files of all coverage maps are on http://hermes.sussex.ac.uk/ as will be new coverage maps as the data 
are taken. 


scan ‘legs’) is 348", nearly the full projected array size (see 
Figure 3.2 of the SPIRE Observers’ manual, V2.4). SPIRE 
observations use two near-orthogonal default scan angles i.e. 
±42.4°. 

Multiple map repeats were required to integrate down 
to the flux limit in each level. These repeats were performed 
with as much cross-linking as possible (i.e. with similar num- 
bers of scans in quasi-orthogonal directions), to enable map- 
ping with the presence of low-frequency drifts and redun- 


dancy for the removal of any problematic scans. We used 
the nominal SPIRE scan rate of SC/'s^ 1 for these fields. 

Where long observations had to be split we aimed to 
cover the whole field on separate occasions (rather than di- 
viding the field and subsequently building a mosaic) to give 
redundancy and maximal cross-linking. 

The Lockman SWIRE and CDFS SWIRE observations 
in Level 5 (#27 and 28) were motivated by the study of 
extragalactic background fluctuations. 

These observations required the rapid scanning using 







8 S.J. Oliver et al. 


the SPIRE fast scan rate at GCf's -1 to minimize the effects of 
low-frequency drifts and increase redundancy. The scanning 
angles and scan leg offsets are the same as for the nominal 
scan rate. 

The knee frequency is that at which the power of the 
correlated fluctuations (primarily from the thermal drifts) 
equates to the white noise. The design goal for the SPIRE 
detectors was for the knee to be at 30 mHz (with a require- 
ment of 100 mHz) but the in-flight performance is much bet- 
ter and by using the thermometer signals to de-correlate the 
drifts knee frequencies of 1-3 mHz can be recovered (Griffin 
et al. 2010). The drift is correlated across the detector array 
(139 bolometers at short wavelengths) and so the effective 
knee frequency for maps is higher. Assuming the knee fre- 
quency to be 30 mHz thermal drift effects would impact on 
a spatial scale of 33' (for the fast scan rate) compared to 17' 
for the nominal scan rate. 

Levels 5 and 6 (#29-41 and 22B) are being mapped 
with the SPIRE-PACS parallel mode. This mode is de- 
scribed in detail in the SPIRE-PACS Parallel Mode Ob- 
servers’ Manual. 4 Parallel mode maps the sky simultane- 
ously with both instruments. The SPIRE detector sampling 
rate is reduced from 18.2 Hz to 10 Hz in this mode, which 
has a negligible impact when scanning in the slow (20”s _1 ) 
mode. The PACS instrument In the blue channel we used 
the PACS Blue2 85-125 /jrn filter (rather than the 60-85) for 
maximum sensitivity. We used the 20 ,, s~ 1 scanning mode as 
the GCf's -1 mode was not suitable for PACS as the beam is 
degraded by up to 30 per cent (Poglitsch et al. 2010, and 
Table 4). 

The parallel mode achieves the combined PACS and 
SPIRE sensitivities more efficiently for large areas than ob- 
servations using each instrument in turn. Scan directions al- 
ternate between nominal and orthogonal for maximal cross 
linking. 

The Level 7, HeLMS, observations (#42) exploited the 
ability of the SPIRE to make long (20 deg) scans at the fast 
(GCf's -1 ) scan rate. These were interleaved in a cross-like 
configuration to give duplicate coverage in a near-orthogonal 
direction. The resulting 270 deg 2 maps are thus optimised 
for studying fluctuations on the largest possible scale. 

All PACS-only observations (Levels 3-4, #20, 21, 25, 
26) were carried out using the scan mapping mode. This 
mode is described in detail in the PACS Observers’ Manual. 5 

The noise of the PACS bolometer/readout system has 
a strong 1// component (Poglitsch et al. 2010) and obser- 
vations need to be modulated on a time-scale of 1-5 Hz. We 
used the 20 ,, s _1 scan rate in which the beam has FWHM 
~6.8" or ^11.3" in the two bands we use (see Table 4), i.e. 
sources are modulated on ~2-3 Hz time-scale. Faster scan 
rates (e.g. GO^s -1 in parallel mode) would have introduced 
significant beam smearing of around 30 per cent (Poglitsch 
et al. 2010, and Table 4). 

We alternated orthogonal scan directions to minimise 


4 The SPIRE-PACS Parallel Mode Observers’ Man- 

ual is available from the Herschel Science Centre 
http:/ /herschel. esac.esa.int/Docs/PMODE/html/paralleLom.html 

5 The PACS Observers Manual is available from the Herschel Sci- 
ence Centre http://herschel.esac.esa.int/Docs/PACS/html/pacs 
om.html 


correlated noise, i.e. correlations arising from asymmetric 
transient detector responses to sky signal. 


2.4 Dithering 

Moving the array on successive scans so that different pixels 
or bolometers trace different parts of the sky (dithering) im- 
proves the quality of the data in a number of ways. Dither 
steps of more than one detector will reduce correlated noise 
arising when the same detector crosses the same patch of 
sky on a short timescale. Dithering on large scales will also 
increase uniformity by distributing dead/noisy pixels across 
the maps. Dithering at sub-detector scales can possibly lead 
to some improvement in resolution if the point spread func- 
tion is not fully sampled and (in the case of SPIRE) further 
reducing the impact of the sparse filling of the focal plane. 

For PACS-only observations we implemented a dither- 
ing pattern. For each scan we requested an offset with re- 
spect to our nominal target position with offsets defined on 
a grid with spacing (0 ,, ) ±7. 5", ±10.5”). This provides sam- 
pling at sub-pixel and sub-array scales. 

For SPIRE we modelled the scan pattern of good de- 
tectors and investigated dithering patterns that reduced the 
variation in sensitivity to point sources (for details see Ap- 
pendix B). We found that for a given number of repeats, 
N, offsetting by a fraction 1/A of the scan leg separation 
between repeats was usually close to optimal. Exceptions to 
this would be cases where the resulting step size coincided 
with the projected bolometer spacing, however, none of our 
patterns resulted in that coincidence. This also provided a 
good de-correlation of the noise. The disadvantage of these 
large dither steps is that the coverage declines at the edges 
of the map. However, for our large maps this is not a ma- 
jor penalty. Since each SPIRE-only observation consisted of 
two scans one at each of the near-orthogonal SPIRE scan 
angles we set an offset in both directions at once. We ar- 
ranged these offset pairs in a square pattern to minimise the 
edge effects. This dithering was not done for observations 
taken during the Science Demonstration Phase, but was im- 
plemented afterwards. The contrast in the coverage maps 
between dithering and not dithering can be seen in Fig. 4. 


2.5 Sensitivity 

To estimate the sensitivity of our survey design we use the 
Herschel Observation Planning Tool, HSPOT v5.1.1. For our 
survey scanning patterns we compute the 5<r instrument sen- 
sitivity (ignoring confusion noise). The HSPOT sensitivities 
are tabulated in Table 3 and their implications for Herschel 
surveys in Table 5. 


2.6 Economies from nesting 

We have designed our survey starting at the widest, shal- 
lower tier and building up the deeper tiers. Thus a small 
field tier nested within a shallower tier needs fewer repeats 
to reach the required depth. This improves the overall survey 
efficiency, because observations of small fields are relatively 
inefficient due to the overheads associated with telescope 
turn-arounds. 
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The current coverage of the nested fields around CDFS 
is illustrated in Fig. 6. 

The nesting of fields is indicated in columns 5 and 6 
in Table 5. and the sensitvities in Table 5 take this into 
account. E.g. UDS-HerMES at Level 3 (#21) includes 12 
PACS scans from UDS Level 4 (#25), in addition to the 
25 from Level 3, giving a total of 37 as well as 14 SPIRE 
nominal scans from UDS Level 4 (#23), four Parallel scans 
from XMM- VIDEO at Level 5 (#32) and two Parallel scans 
from Level 6 XMM-LSS SWIRE (#36). 

2.7 Total Time 

The total time allocated for HerMES is 909.3 hours. This 
comes from the Guaranteed Time awarded to the SPIRE 
instrument team (850 hr) one of the Herschel Mission Sci- 
entists (M. Harwit, 10 hr) and members of the Herschel Sci- 
ence Centre (B. Altieri, L. Conversi, M. Sanchez Portal and 
I. Valtchanov, 40hr). ESA also effectively contributed 9.3 
hours as we agreed for our Abell 2218 observations in Sci- 
ence Demonstration Phase to be made public immediately 
and so were not charged for these observations. 


2.8 Special requirements and constraints 

The Herschel observatory is performing very close to specifi- 
cations and our survey design is very similar to the one pro- 
posed. However, some changes and compromises have been 
made on the basis of post-launch experience. 

Early in the mission there was a constraint that parallel 
mode observations could not exceed 2 15 s, as this exceeded 
the limit of one software counter. Since each parallel mode 
observations was already a single-scan they were as shallow 
as could be done at that scan rate so this required us to 
split some of the Level 5 and 6 fields into smaller fields, 
compromising the uniformity of the data. The impact of 
this on the coverage for the XMM-LSS and Bootes fields is 
shown in Fig. 5. The planned AKARI deep field south (#41) 
and ELAIS SI (#39) fields required only slightly more time 
than 2 15 s, and so we chose to reduce the field size rather 
than split the field. 

Where the orientation of the SPIRE data with respect 
to complementary data was particularly important we con- 
strained the observations to align with them. Solar avoid- 
ance constraints meant that it was not possible to align the 
SPIRE observations of XMM-LSS SWIRE (# 36) and COS- 
MOS optimally with the Spitzer data and PEP data, re- 
spectively. For XMM-LSS SWIRE we observed a larger field 
containing the Spitzer data, while for COSMOS we observed 
a larger shallower field, COSMOS HerMES (#22B), con- 
taining the planned PEP PACS observations and a smaller 
deeper field (COSMOS, #22), which does not fully cover the 
PACS observations. 

The Spitzer - SERVS and VISTA- VIDEO surveys were 
approved after HerMES and designed with reference to Her- 
MES. So, almost all the SERVS and VIDEO fields were in- 
cluded in our Level 5 observations. However, the SERVS and 
VIDEO field in ELAIS SI was not quite within our planned 
observations, which were only at Level 6. We thus included 
additional deeper observations covering the SERVS/ VIDEO 
field (#39B). 



Figure 6. Map of square root of number of effective number 
of bolometers samples per pixel for SPIRE 250 /im blank field 
observations of the CDFS region, which includes our GOODS- 
S, ECDFS and CDFS-SWIRE observations (#13,15,27,33). The 
parallel mode samples (#33) have been scaled by the relative 
sampling rates, 18.6/10, to give the effective number of samples 
they would have had if the observation had been carried out with 
SPIRE large-map mode with the same exposure time per pixel. 
A region of uniform coverage for each of the independent sets 
is marked with a rectangle. N.B. the total coverage drops off in 
the north-eastern corner of the largest rectangle (delimiting #33) 
due to the coverage coming from the boundaries of the large-map 
mode observations (#27) but is uniform in a coverage map built 
from #33 data alone. 


Our initial Science Demonstration Phase (SDP) obser- 
vations of Abell 2218 used ‘Large Map’ mode as this was 
the only way of doing scan mapping. We changed our deep 
cluster observations to the ‘Small Map’ mode once the mode 
was available. 

The P(D) results of Glenn et al. (2010) successfully 
probed the number counts well below the confusion limit, 
reducing the motivation for exceptionally deep cluster ob- 
servations. We have thus reduced the number of repeats. 

Due to an error in entering the AOR one parallel obser- 
vation scan of ELAIS SI SWIRE (#39) was accidentally ob- 
served with the shorter wavelength 60-85 /im channel rather 
than the 85-125 pm channel. 

The PACS sensitivity of lOmJy (5-c in 1 hr) in the 
85-125 channel was significantly less than the pre-launch 
estimate (5.3 mjy, PACS Observers’ manual vl.l) and we 
removed our planned PACS observations of the VVDS field 
(# 26). 

To extend the fluctuation science goals and increase the 
Herschel discovery space for rare objects including gravita- 
tionally lensed systems, we added the HeRMES Large-Mode 
Survey (HeLMS), a wide, SPIRE only, tier of 270 deg 2 tak- 
ing around 100 hours. This exploits the ability of SPIRE to 
cover wide areas close to the confusion limit. This additional 
level is indicated in Table 5. 
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Figure 5. Maps of the number of bolometer samples per pixel in SPIRE 250 fun blank field observations from Level 6. From the left they 
are XMM-LSS-SWIRE (#39), Bootes NWDFS (#37) taken early with conservative overlap) and FLS (#40, from SDP). All are parallel 
mode observations with a nominal coverage of two scans. Overlaps produce a maximum coverage of four scans in XMM-LSS-SWIRE 
and eight in Bootes. 


5-cr sensitivities [mjy \/N SC an 
Mode Rate Step at Wavelength [/im] 



["s- 1 ] 

["] 

100 

160 

250 

350 

500 

SPIRE 

30 

348 



64 

53 

76 

SPIRE 

60 

348 



91 

75 

108 

PACS 

20 

55 

42 

80 




Parallel 

20 

168/155 

71 

135 

37 

30 

44 

Parallel 

60 

168/155 

122 

232 

63 

53 

75 


Table 3. Point source sensitivities for different Herschel observ- 
ing modes. Scan rates are given for each mode, we also tabulate 
the step size between successive scan legs (pre-determined for 
SPIRE and parallel mode but user-defined for PACS). In parallel 
mode the step size are different for maps built by scanning in 
each of the two “orthogonal” directions. 5-cr sensitivities in units 
of [mjy \/N sc an] for a single scan are estimated from the HSPOT 
v5.1.1. Modes below the line are not used by HerMES but by 
other Key Program surveys. 

Beam FWHM ["] 


Mode 

Rate 

["s- 1 ] 

100 

at Wavelength [/im] 
160 250 

350 

500 

SPIRE 

PACS 

30/60 

20 

6.8 

11.4 

18.2 

24.9 

36.3 

Parallel 

20 

6.8 

11.4 

18.2 

24.9 

36.3 

Parallel 

60 

7.0x12.7 

11.6x15.7 

18.2 

24.9 

36.3 


Table 4. Beam sizes for different Herschel observing modes. Scan 
rates are given for each mode. The FWHM of the beams in 
units of ["] are taken from SPIRE and PACS Observers’ Man- 
uals V2.4/V2.3 (respectively). Where two values are given these 
are the major and minor axes, when the ellipticity is less than 
15 per cent the geometric mean of the two is quoted. The SPIRE 
beam is not known to vary significantly with scan rate and only 
one value is given. Modes below the line are not used by HerMES 
but by other Key Program surveys. 

2.9 Observations 

Our first observation was carried out on 12th September 
2009. This was the first half of our SPIRE observations of 
Abell 2218 (#1) and the resulting map from all the data 
is shown in Fig. 7. This was part of the Herschel Science 
Demonstration Phase (SDP). Our SDP observations were 


designed to exercise most of the modes that were to be used 
in the full survey, and the SPIRE observations are described 
in Oliver et al. (2010b). This includes the observations of 
GOODS-N (#14) (Fig. 8). The SDP observations concluded 
on 25th October 2009; AORs are available under the pro- 
posal ID SDP_soliver_3. 

The program is now being carried out 

as part of the Routine Phase (proposal ID 

KPGT_soliver_l) and is expected to be completed 
during 2011. The current ESA schedule is on 

herschel . esac . esa. int/observing/ScheduleReport .html 
and the observing log can be followed on 

herschel . esac . esa. int/observing/LogReport .html. 

2.10 Comparison with other Herschel Surveys 

HerMES was planned alongside the ‘PACS evolutionary 
Probe (PEP)’ survey (Proposal ID KPGT_dlutz_l, e.g. Lutz 
et al. 2011). Since then there have been a number of related 
Key Project surveys carried out in Open Time. There have 
also been a few surveys being undertaken in Open Time but 
not as Key Projects. These programmes are listed in Table 6 

The cumulative area of all major Herschel-SPIRE extra- 
galactic Key Program surveys as a function of instrumental 
noise (taken from Table 5) and for the HerMES fields is 
shown in Fig. 9. 

It is striking to compare the Hersc/iel-SPIRE sub- 
millimetre surveys with previous sub-millimetre surveys. To 
do this we have explored the sensitivity of surveys to a 
canonical galaxy with a modified blackbody spectral en- 
ergy distribution with emissivity, /3 = 1.5, and temperature 
T = 35 K. These are shown in Fig. 10. 

3 EARLY AND ANTICIPATED SCIENCE 
3.1 Confusion limits 

An important consideration in design of HerMES was the 
impact of source confusion at SPIRE wavelengths, i.e. the 
limited ability to separate individual sources due to the res- 
olution of the telescope and the number density of sources. 
We define confusion noise to be the standard deviation of 
the intrinsic variations in a map on the scale of the beam 
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Fields 

Area 
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[deg 2 ] 

Cummulative 

Observations 

PACS 

SPIRE 

5 -a noise level (for band in 
110 160 250 350 

[mjy] 

/cm) 

500 

Abell 2218 

0.0050 
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0.1 

p 

1 

4.1 

7.9 

6.4 

5.3 
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Abell 1689 
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0.11 

p 

2 

3.6 
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9.2 

7.7 

11.0 

8 Targets 

0.04 

0.04 

0.15 

p 
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10.9 

9.2 

7.7 

11.0 

2 Targets 

0.03 
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0.18 

p 
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13.9 

11.6 

16.7 
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0.18 

0.18 

0.36 

E 

E 

6.1 

11.7 

14.2 

11.9 

17.1 

GOODS-N 

0.042 

0.042 

0.04 

G,P 

G,14 

2.2 

4.1 

3.8 

3.1 

4.5 

GOODS-S 

0.11 

0.087 

0.13 

G,P,33 

13,15,27,33 

2.1 

2.9 

4.3 

3.6 

5.2 

GOODS-S 

0.012 

0.012 

0.14 

G,P,33 
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1.1 

2.1 

4.6 

3.8 

5.5 

GOODS-S 

0.018 

0.0060 

0.15 

G,P,33 

13,15,27,33 

1.6 

3.0 

4.6 

3.8 

5.5 
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0.023 

0.0060 

0.15 

G,P,33 
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2.0 

2.0 
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P 
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8.0 

6.6 
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ECDFS 

0.25 

0.14 
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14.5 

8.0 

6.6 
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GOODS-N 

0.25 

0.208 

2.5 

P 

14 
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8.9 

8.2 

6.8 

9.9 

Lockman-East 

0.25 

0.25 

2.75 

P 

18,18B,28B,34,28 

6.5 

12.3 

9.6 

7.9 

11.5 

Lockman-North 

0.25 

0.25 

3.0 

20,20B,34 

19,28B,34,28 

7.4 

14.1 

10.6 

00 

bo 

12.7 

Groth Strip 

0.25 

0.25 

3.25 

P,29 

17,29 

7.1 

13.6 

10.7 

8.9 

12.8 

UDS HerMES 

0.25 

0.25 

3.5 

21,25,32,36 

23,25,32,36 

6.8 

12.9 

11.2 

9.3 

13.4 

UDS 

0.7 

0.7 

4.2 

25,32,36 

25,32,36 

11.2 

21.4 

11.2 

9.3 

13.4 

VVDS 

2.0 

2.0 

6.2 

25,320,36 

25,320,36 

28.8 

54.9 

11.2 

9.3 

13.4 

CDFS SWIRE 

11.4 

11.1 

17.3 

33 

27,33 

31.5 

60.2 

12.7 

10.5 

15.2 

Lockman SWIRE 

16.1 

15.6 

32.9 

34 

28,28B 

35.3 

67.3 

13.6 

11.2 

16.2 

EGS HerMES 

2.7 

2.5 

35.4 

29 

29 

26.6 

50.8 

13.8 

11.3 

16.4 

Bootes HerMES 

3.3 

3.3 

38.6 

30,37 

30,37 

26.6 

50.8 

13.8 

11.3 

16.4 

ELAIS N1 HerMES 

3.3 

3.3 

41.9 

31,35 

31,35 

26.6 

50.8 

13.8 

11.3 

16.4 

ELAIS SI VIDEO 

3.7 

3.7 

45.6 

39B,39 

39B,39 

28.8 

54.9 

14.9 

12.2 

17.8 

XMM-LSS VIDEO 

7.7 

5.0 

50.6 

32,32B,32C,36 

32, 32B, 32C,36 

28.8 

54.9 

14.9 

12.2 

17.8 

COSMOS Hermes 

4.4 

2.4 

53.0 


22B 



15.9 

13.3 

19.1 

ELAIS N2 SWIRE 

7.9 

7.9 

60.9 

41 

41 

49.9 

95.1 

25.8 

21.2 

30.8 

FLS 

6.7 

6.7 

67.6 

40 

40 

49.9 

95.1 

25.8 

21.2 

30.8 

ADFS 

7.5 

7.5 

75.1 

38 

38 

49.9 

95.1 

25.8 

21.2 

30.8 

ELAIS SI SWIRE 

7.9 

4.2 

79.2 

39 

39 

49.9 

95.1 

25.8 

21.2 

30.8 

ELAIS N1 SWIRE 

12.3 

9.1 

88.3 

35 

35 

49.9 

95.1 

25.8 

21.2 

30.8 

Bootes NDWFS 

10.6 

7.3 

95.6 

37 

37 

49.9 

95.1 

25.8 

21.2 

30.8 

XMM-LSS SWIRE 

18.9 

15.0 

110.6 

36 

36 

49.9 

95.1 

25.8 

21.2 

30.8 

Various 

570.0 

570.0 

681.0 

A 

A 

86.3 

164.0 

44.5 

37.1 

53.0 

SPT 

100.0 

100.0 

781.0 


S 



45.3 

37.5 

54.1 

HeLMS 

270.0 

270.0 

1051.0 


42 



64.0 

53.0 

76.5 


Table 5. HerMES survey with sensitivities in the context of other survey programmes being undertaken by Herschel. The “observations” 
columns refer to the AOR set numbers of Table 1 for HerMES or for other Key Programmes we use: “E” for Egami cluster programme, 
“G” for GOODS-H, “P” for PEP, “A” for H- ATLAS and “S” for SPT (see Table 6). The sensitivities are estimated consistently using 
HSPOT v5.1.f. These are single pixel sensitivities and ignore the benefits of matched filters, particularly for unconfused fields, e.g. 
H-ATLAS quote empirical 5-cr sensitivities of 105, 140, 32, 36, 45 mjy for the five wavelengths so the sensitivities in this Table should be 
scaled by 1.22, 0.85, 0.72, 0.97, 0.85 to obtain a consistent comparison with H-ATLAS. The sensitivity of HerMES observations have been 
calculated including data from shallower tiers as described in the text. Other surveys are treated independently. Cluster observations are 
listed before blank fields. The fields are ordered in increasing 250 (im flux limit then right ascension. The area is defined by the PACS 
observations for Levels 1-4 (above the second horizontal line), otherwise we use iL 00(j from Table 1 or f2 n0 m for HeLMS. We tabulate 
three areas: the nominal area for each field; the ‘doughnut’ area which excludes any deeper sub-fields within; and the cumulative area 
of all fields higher in the table. The 5-<r confusion noise (after 5c cut) from Nguyen et al. (2010) is 24.0, 27.5, 30.5 mjy (at 250, 350 
and 500 ttm), approximately the Level 6 depth. GOODS-S also has PACS data not listed here at 70//rn over 0.11 deg 2 to a 5-cr depth of 
1.9 mjy. 


due to all point sources . We planned our survey with refer- 
ence to several number count models (Lagache et al. 2003; 
Le Borgne et al. 2009; Franceschini et al. 2010; Pearson & 
Khan 2009; Xu et al. 2003). We used these models to esti- 
mate the fluctuations in a map which at the 4 -a level were 
1.6 ± 0.9, 10.6 ±3.1, 26.3 ± 6.3, 32.5 ± 7.5 and 30.0±7.5 mjy 
at 100, 160, 250, 350 and 500 /im respectively. The uncer- 
tainties come from the scatter between models. The SPIRE 


confusion noise estimates compare very favourably with the 
fluctuations in our maps as calculated by Nguyen et al. 
(2010) with 5cr = 24.0,27.5,30.5 mjy at 250, 350 and 
500 fi m, respectively after cutting maps at 5cr. This is per- 
haps fortuitous given that the model counts do not fit the 
observed counts very well in detail (e.g. Oliver et al. 2010b; 
Glenn et al. 2010) but may be because the models had been 
constrained to fit the infrared background. 
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Figure 7. Three colour Herschel-SPIHE image of the central 4/ x 4' of the galaxy cluster Abell 2218. The left-most panels show the 
single band images of the cluster, while the central panel shows a three colour image generated by resampling the single band images and 
their flux scalings to a common pixelization. The centre of the cluster is marked with the cross hairs and a V bar is shown for scaling; 
north is toward the top of the page. The orange object to the south-east and white object to the south-west of the cluster are images 
of the multiply imaged sub-mm source studied in detail by e.g. Kneib et al. (2004); this source has been identified to lie at z = 2.516 
though due to the complex mass structure of this cluster each image is magnified by a different amount. In the SPIRE bands this source’s 
integrated flux densities are measured to be {170, 197, 231} mJy, corresponding to background flux densities of {11.7, 13.5, 15.4} mJy. 
The varying colour of the images suggests that different regions of the source galaxy are being imaged to different points in the map. 
In addition, the known z = 4.04 sub-mm source is seen as the pink object just to the east of the cross hairs (Knudsen et al. 2009). The 
other objects scattered through the image are more typical z 1 sources with SEDs peaking shortward of 250 /xm. 


Call 

Title 

Proposal ID 

Time 

[hr] 

Reference 

Key GT 

Herschel Extragalactic Multi-tiered Survey (HerMES) 

KPGT_soliver_l 

806 

this paper 

Key GT 

PACS evolutionary Probe (PEP) 

KPGT_dlutz_l 

655 

Lutz et al. 2011 

Key OT 

The Cluster Lensing Survey 

KPOT_eegami_l 

292 

Egami et al. 2010 

Key OT 

The Herschel Astrophysical Terahertz Large Area Survey’ (H- ATLAS) 

KPOT_seales01_2 

600 

Eales et al. 2010a 

Key OT 

The Great Observatories Origins Deep Survey (H-GOODS) 

KPOT_delbaz_l 

363 

Elbaz et al. 

OT1 

The Herschel - AKARI NEP Deep Survey 

OTl_sserje01_l 

74 

Serjeant et al. 

OT1 

A deep PACS survey of AKARI-Deep field south’ 

OTl_ttakagi_l 

35 

Takagi et al. 

OT1 

SPIRE Snapshot Survey of Massive Galaxy Clusters 

OTl.eegami 

27 

Egami et al. 2010 

OT1 

Measuring the Epoch of Reionization 

OTl_jcarls01_3 

79 

Carlstrom et al. 

GT2 

HerMES Large Mode Survey 

GT2_mviero_l 

103 

Viero et al. & this paper 


Table 6. Herschel blank field and cluster lens surveys carried out as Key Programmes or ordinary programmes under Guarenteed Time 
(GT) or Open Time (OT). 


HerMES 13 


250|jm 

9 - : 

••iff' * / 9 j • # 

GOODS-N- 

350pm 

^ a v ft * ^ _ ft 1 » 

_ * 0 * ji % 1 

• j 1 v ' 4 f • 

— .% • - a ^ # » 

. 1 1 • V « ft. t ft • ft <i 

•» ft.i >• „ • »,.f ^ <2 ; * 

• \ •« 

r V 

• 1 ft § • 'I m 

r • k f . # • • ft _ w 

* ft % m • 

* • . • 

-ft • * 

ft _ ft# ^ ^ M A ^ ft ftA 

500pm 

■ ■ » ^ ^ W - ft ft . ^ 

ft 

ft • • • ft! 

» ^ft % Y ft • * ^ A ft ^ \ 

|r A ^ • r ft" 

<• * 

* • 


• ' \ ft- • J * 

'•ft.-MC 7 1 >i . • 7* 

• ^ ft . • 

lOarcmin * 

• ◄ ► • 


Figure 8. Three colour Herschel-SPIRE image of the GOODS-North region. This is a sub-set of our GOODS-N observation. The 
left-most panels show the single band images of the cluster, while the central panel shows the three colour image. 
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Figure 9. Cumulative area against 5-cr instrumental noise level 
at 250 for the HerMES blank-field surveys with SPIRE. The 
colour-coding breaks this down into individual survey fields. 


We planned for the survey to have a substantial area 
(providing SDSS-like volumes) at the confusion limit, but 
with some regions well below the confusion limit in very well 
studied fields, to exploit techniques for mitigating confusion 
using high signal-to-noise data. 
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Figure 10. Luminosity limit verses redshift for submm surveys 
to date. The luminosity limit was calculated assuming a modified 
blackbody of 35K at z = 2. (References for the points are as 
follows: SCUBA - Hughes et al. 1998; Scott et al. 2002; Coppin 
et al. 2006, MAMBO - Greve et al. 2004; Bertoldi et al. 2007, 
BOLOCAM - Laurent et al. 2005, AzTEC - Perera et al. 2008; 
Austermann et al. 2010; Scott et al. 2010b; Aretxaga et al. 2011, 
LABOCA - Weifi et al. 2009, SPT - Vieira et al. 2010; Williamson 
et al. 2011, BLAST - Devlin et al. 2009, SPT SPIRE - Carlstrom 
et al.) 
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Figure 11. The comoving infrared luminosity density, a proxy 
for the star formation history of the Universe, from Seymour 
et al. (2010). In grey shading is the Spitzer view (following Le 
Floc’h et al. 2005) showing the contribution of low IR luminosity 
galaxies (solid), luminous IR galaxies (LIRGs, cross-hatched), and 
rapidly-evolving ultra-luminous IR galaxies (ULIRGs, hatched), 
to the total co-moving IR energy density. Different symbols give 
subdivision by dust temperature within each luminosity class. 

3.2 Science above the confusion limit 

3.2.1 Direct determination of the total far infrared 
luminosity function and its evolution 

Our primary goal has been to determine the total far in- 
frared luminosity function and subsequently the bolometric 
luminosity of galaxies over the redshift range 0 < 2 < 3. 
For this analysis we use galaxies detected in Herschel im- 
ages combined with extensive multi-wavelength data to de- 
termine photo-as where spectroscopic redshifts are not yet 
available. 

Our first results on exploration of the full far infrared 
SED are given by Elbaz et al. (2010); Rowan-Robinson et al. 
(2010); Hwang et al. (2010) and Chapman et al. (2010). El- 
baz et al. (2010) combined photometry from PACS (from the 
PEP program) and SPIRE (from HerMES). We found that 
the total far infrared luminosity estimated from extrapola- 
tions of Spitzer 24 pm data agreed well with direct measure- 
ments from Herschel at lower redshift but underestimated 
the power at higher redshifts (as also seen by Nordon et al. 
2010). In that work the longer wavelength (SPIRE) band 
measurements departed from the model SEDs at lower red- 
shift. This was explored further by (Rowan-Robinson et al. 
2010), showing that the SPIRE results for some galaxies 
could be explained with a cold dust component not nor- 
mally included in canonical templates. Indeed, when simply 
characterising the SEDs by their effective dust temperature 
we have shown that the SPIRE detected galaxies cover a 
broad range of temperatures (Hwang et al. 2010; Magdis 
et al. 2010) and thus capture warm objects like the ‘Op- 
tically Faint Radio Galaxies’ missed by ground-based sub- 
millimetre surveys (Chapman et al. 2010). 

We have already determined our first measurements of 
the local luminosity functions at 250, 350 and 500 pm to- 
gether with a total infrared (8-1000 pm) function, finding a 
local luminosity density of I.SIq'j x 10 8 Lq Mpc -3 (Vaccari 



Figure 12. A slice of the dark matter in the Millennium Simula- 
tion of the Universe, seen today (Springel et al. 2006). Overlayed 
are the footprints of some of our fields, showing how much of this 
slice they would sample at z = l. This thin slice exaggerates the 
effect but illustrates that to overcome sampling variance and to 
probe a full range of environments we need multiple, large fields. 


et al. 2010) and showing that the 250 pm function evolves 
strongly to z ~ 1 (Eales et al. 2010b), similarly to earlier 
studies at shorter wavelengths. Future analysis (in prepara- 
tion) will study wider areas with more and better ancillary 
data and extend these results to higher luminosities, higher 
redshifts and model the relative contribution of AGN and 
star-formation to the bolometric emission, as well as ex- 
ploring the relation between the infrared luminosities and 
the stellar properties probed at optical, NIR and UV wave- 
lengths. 


3.2.2 Star- formation and environment 

Environment on various scales plays an important role in 
the process of galaxy formation. Perhaps the most striking 
observational evidence is that clusters today have a much 
higher fraction of early-type galaxies than is found in the 
field. Likewise the successful physical models of galaxy for- 
mation predict a very strong co-evolution between galaxies 
and dark-matter halos. 

There are many ways of determining the role of environ- 
ment observationally: one can directly examine the galaxy 
properties (e.g. the SFR distribution functions) in different 
environments; one can explore the environments of galaxies 
in different luminosity classes; one can use the clustering of 
particular galaxy populations to infer the mass of the dark 
matter halos in which they are located, to relate these to 
their present-day descendants; or one can directly use the 
structure in the maps to constrain such models. All these 
methods have the same basic requirement, a volume suffi- 
ciently large to sample enough of the environments of in- 
terest, and sufficiently deep to constrain the populations of 
interest. A simulation in Fig. 13 shows that we could dis- 
criminate different halo mass hosts for different sub-classes 
of galaxies and compare the clustering of the FIR galaxies 
with quasars from optical studies. 

First results on the clustering of HerMES galaxies were 
given by Cooray et al. (2010), indicating that the HerMES 
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Figure 13. Evolution of co-moving correlation length, 7*0 , with 
redshift. Solid lines show the predicted clustering amplitude of 
haloes of given mass. We have simulated data for the clustering of 
LIRGs (red), ULIRGs (green) and HLIRGs (magenta), assuming 
they inhabit halos of mass 10 12 , 10 13 and 1O 135 M0 respectively. 
The simulation is for our 250 pm surveys at Level 5 (square) and 
Level 6 (triangle) . For comparison we show quasar clustering from 
Croom et al. (2005) as stars and SCUBA galaxies from Blain et al. 
(2004) as orange circles. Spitzer sources from Farrah et al. (2006) 
are shown as blue circles and blue shaded regions extrapolate 
those objects to their progenitors and descendants. 


sources with S 250 > 30 mjy (at z ^ 2) were in dark matter 
halos with masses above (5 ± 4) x 10 12 Mq. 

Clustering can also be used in other ways. A recent 
cross-correlation analysis indicates that there is a correlation 
between HerMES sources at s ~ 2 and foreground galaxies 
from SDSS at 2 ~ 0.2 and SWIRE at z ~ 0.4 (Wang et al. 
2011). While some of this signal can be attributed to the 
intrinsic correlation of galaxies in the overlapping tails of 
the redshift distributions, there is clear evidence for a signal 
arising from the amplification of the HerMES source fluxes 
by lensing from foreground galaxies. 


3.2.3 Extreme galaxies 

Rare objects provide challenges for theories and may ex- 
pose important but transitory phases in the life-cycle of 
galaxies. The very wide surveys, in particular, will discover 
many exotic objects, which are prime targets for ALMA. 
Galaxies with extremely high star formation rates would 
be difficult to explain with some models of galaxy forma- 
tion. Limited area sub-millimetre surveys have already dis- 
covered small samples of galaxies with very high star for- 
mation rates > 1000MQyr _1 e.g. SMM J02399-0136 (Ivison 
et al. 1998), GN20 (Daddi et al. 2009; Borys et al. 2003) 
and MIPS J142824. 0+352619 (Borys et al. 2006). By map- 
ping large areas at the wavelengths where re-emission from 
star formation peaks, we will be able to quantify the number 
density of systems of ~ lOOOMQyW 1 and determine whether 
there are any systems with even higher star formation rates. 
Even individual examples of such systems would be impor- 
tant as extreme astrophysical laboratories and would pro- 
vide fruitful targets for new facilities, especially ALMA. 

A primary search tool will be the SPIRE colours. 
Searches have already revealed many galaxies (Schulz et al. 


2010) with very red colours S 25 o/<S 350 < 0.8 and with flux 
densities above 50 mjy. These may be a mix of intrinsically 
cool galaxies at lower redshift, and galaxies at high redshift, 
including some that are lensed by foreground galaxies. 

3.2.4 Lensed Systems 

Lensed systems are interesting because, although lensing is 
a rare phenomenon, they provide a magnified view of more 
common, relatively normal, but distant galaxies, which can 
then be easily studied. An example of a lensed source found 
in early HerMES data is HERMES J105751. 1+573027, a 
2 = 2.957 galaxy multiply lensed by a foreground group of 
galaxies. Coupled with a lensing model derived from high- 
resolution observations (Gavazzi et al. 2011), the magnifi- 
cation and large image separation allowed us to investigate 
the continuum SED from the optical to far-IR (Conley et al. 

2011) , as well as model the CO line excitation (Scott et al. 
2011) and study the gas dynamics (Riechers et al. 2011). 

3.3 Science below the confusion limit 

The deepest observations at SPIRE wavelengths suffer sub- 
stantial confusion noise due to faint unresolved galaxies, and 
are limited in their ability to define true luminosities, SEDs 
and physical conditions within the most active galaxies dur- 
ing the peak epoch of galaxy formation at redshift z ~ 2. 
We will investigate and employ super-resolution techniques, 
e.g. CLEAN (Hogbom 1974) or matched filtering (Chapin 
et al. 2011). However, as argued in Oliver (2001), we expect 
the gains from blind image deconvolution techniques to be 
modest except at the very highest signal-to-noise ratios. 

One approach to combat the problem is to study iso- 
lated sources as we have discussed in Elbaz et al. (2010); 
Brisbin et al. (2010) and Schulz et al. (2010), however, we 
are pursuing many other mitigating techniques. 

3.3.1 Ultra-deep far-infrared galaxy surveys from imaging 
of rich clusters of galaxies 

Rich clusters can be used as tools to mitigate this effect, al- 
lowing high-redshift galaxy formation to be investigated by 
the gravitational magnification of the primordial galaxies 
behind the cluster. This has been demonstrated at relevant 
wavelengths by Smail et al. (2002), Cowie et al. (2002), Met- 
calfe et al. (2003), Chary et al. (2005) and Swinbank et al. 
( 2010 ). 

Gravitational lensing brightens and separates the im- 
ages of all background galaxies within 1-2' of the core of 
the cluster (e.g. Kneib et al. 2004), making individual back- 
ground galaxies easier to detect. This also allows the sources 
of up to about 50 per cent of the otherwise confused and un- 
resolved background radiation to be identified with specific 
galaxies. 

The selected clusters have some of the best archival data 
available, including deep HST ACS/NICMOS images, ultra- 
deep /rjy radio imaging, deep mid-IR imaging from Spitzer, 
and X-ray images from Chandra/ XMM- Newton. The mass 
and magnification profiles are known accurately, from ex- 
tensive spectroscopy of multiply-lensed images (Kneib et al. 
1993). 
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Our observations of 10 clusters will provide about 180 
sources that will allow us to quantify the space density of 
the faintest Herschel galaxies with 10 per cent accuracy. Two 
clusters (Abell 2218 and Abell 1689) were believed, in ad- 
vance, to be relatively free of bright lensed galaxies. This 
was intentional as these were originally intended for very 
deep observations in order to detect of order 10 even fainter 
lensed sources (to determine the counts of Herschel galaxies 
at the 5-mJy detection level, reaching well below the blank- 
field confusion limit) so we wanted to avoid confusion from 
known lensed galaxies. Following modification to our pro- 
gram in the light of analysis of in-flight data only Abell 2218 
was observed deeper than the others. The results from our 
SDP cluster observation of Abell 2218 clearly demonstrate 
that we can detect high redshift lensed galaxies, see Fig. 7. 


3.3.2 Multi-colour 1 -point fluctuation analysis below the 
confusion limit 

Analysis of the fluctuations in the cosmic IR background 
radiation provides unique information on sources too faint 
to be detected individually (see, e.g. Maloney et al. 2005b; 
Patanchon et al. 2009). Our Level 2 and Level 3 fields allow 
us to analyze the fluctuation distribution down to flux densi- 
ties of 2-3 mjy, where much of the background was expected 
to be resolved. By analyzing the fluctuations in all three 
SPIRE wavebands, we can obtain statistical information on 
SEDs. This multi-colour P(D) analysis provides a powerful 
method for distinguishing different number count models, 
thereby constraining the redshifts and emission properties 
of the source population (Fig. 14). This requires very precise 
characterization of the instrument noise for optimal analy- 
sis. 

We undertook a mono-chromatic fluctuation analysis 
using three fields from our SDP data. With that analysis 
(Glenn et al. 2010) we reached a depth of 2 mjybeam -1 , 
significantly deeper than any previous analysis at these 
wavelengths. Modelling this distribution with parameterised 
number counts confirmed the results from resolved sources 
(Oliver et al. 2010b) and was in disagreement with previous 
models. The fits accounted for 64, 60, and 43 per cent of 
the far-infrared background at 250, 350 and 500 pm, respec- 
tively. 


3.3.3 Average SEDs of galaxies contributing to the 
infrared background 

Prior information from shorter wavelength (e.g., 24 pm with 
MIPS) can be used to infer the statistical properties (such 
as source density or SEDs) at longer wavelengths. A more 
promising route to achieving super-resolution results is to 
use prior information on the positions of sources from higher 
resolution data at other wavelengths. This has been demon- 
strated with HerMES data in Roseboom et al. (2010) acheiv- 
ing robust results for source fluxes down to £250 ~ lOmJy. 

A related technique is ‘stacking’, which averages the sig- 
nal from many similar prior sources. In the absence of signif- 
icant correlations the confusion variance would then reduce 
in proportion to the number of prior sources in the ‘stack’. 
Stacking has been successfully applied to Spitzer MIPS data; 
Dole et al. (2006) stacked more than 19,000 24 pm galaxies 
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Figure 14. Simulation of a two dimensional P(D ) analysis, show- 
ing discrimination between models. The x and y axes show the 
pixel intensities (in mjy beam -1 ) in the 250 pm and 350 pm 
bands, respectively. The contours show the number of pixels with 
those intensities, logarithmically spaced. The top panel is for the 
number count model of Valiante et al. (2009), the bottom is for 
the mock catalogues of Fernandez-Conde et al. (2008) based on 
the models of Lagache et al. (2003). The simulations are around 
ten deg 2 and with 1 mjy of Gaussian noise in each band. 



to find the contributions of the mid-IR galaxies to the far- 
IR background (70 and 160 pm). With this technique, they 
gained up to one order of magnitude in depth in the far-IR. 
It appears that a large fraction of the 24 pm sources can be 
statistically detected at longer wavelengths (e.g. Marsden 
et al. 2009). Such an analysis applied to Herschel will allow 
us to extend galaxy SEDs to the FIR/sub-mm to quantify 
the contribution of different populations to the background 
(e.g. Dye et al. 2007; Wang et al. 2006), or to explore the 
star-formation properties as a function of redshift and stel- 
lar mass (e.g. Oliver et al. 2010a). Such procedures might 
use Spitzer 24 pm catalogues and/or the PACS catalogue. 
This type of analysis is critically dependent on the quality 
and depth of the ancillary data, and further motivates our 
choice of very well studied extra-galactic fields. An example 
of this approach is shown in Fig. 15. 

Stacking has already been used in some of our analysis 
(e.g. Ivison et al. 2010; Rigopoulou et al. 2010) and our first 
results analysing the contribution of various prior popula- 
tions to the background through stacking will be presented 
by Vieira et al. (2011). 


3.3.4 Extragalactic Correlations Fluctuations 

A comprehensive fluctuations analysis is an essential com- 
plement to the aspects of our survey allowing us to inves- 
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Figure 15. Example of about 4000 Spitzer IRAC selected ‘Bump- 
s’ sources (i.e., objects with peak emission at 5.8 pm) stacked in 
HerMES SPIRE maps at 250 pm with 6" pixels. The clear detec- 
tion allows one to derive aggregate SEDs of this galaxy popula- 
tion, expected to lie at 2.2 < z < 2.8. 



Figure 16. The angular power spectrum of unresolved 
anisotropies at 350 pm. We model the power spectrum under the 
halo approach and describe non-linear clustering with a halo oc- 
cupation number /3, as shown by the orange curves. We show 
simulated binned errors with SPIRE in the 11 deg 2 Lockman 
Hole L5 field, including both instrument noise and sample vari- 
ance, and removing shot noise from galaxies below the detec- 
tion limit (dashed black curve). For reference, the long-dashed 
and solid blue lines show the noise per multipole for Planck and 
SPIRE, respectively. The green line is the foreground dust spec- 
trum, determined for the same field using dust maps. In red we 
show the residual foreground spectrum after cleaning with multi- 
wavelength data. Even if not removed, dust does not contaminate 
small angular scales, where SPIRE excels. 


tigate the majority population of objects, those below the 
Herschel confusion limit. Using the two shallowest tiers of 
the survey, we can specifically target non-linear clustering 
on angular scales < 10', virtually inaccessible to Planck , 
and where SPIRE is not susceptible to low frequency drifts. 
The clustering of undetected sources produces fluctuations 
on larger spatial scales (Amblard & Cooray 2007; Haiman 
& Knox 2000; Knox et al. 2001) which are expected to be 
brighter (Scott & White 1999) than Poisson fluctuations on 
spatial scales > 1 On large angular scales, background 
fluctuations measure the linear clustering bias of infrared 
galaxies in dark matter halos. On small angular scales, fluc- 
tuations measure the non-linear clustering within individ- 
ual dark matter halos, and the physics governing how FIR 
galaxies form within a halo as captured by the occupation 
number of FIR sources. This halo approach (e.g. Cooray & 
Sheth 2002) will allow us to compare the results of a Herschel 
fluctuations survey with studies at other wavelengths, to ob- 
tain a consistent picture of galaxy clustering and evolution. 
Finally, this fluctuation survey is designed to complement 
surveys by Planck on larger angular scales. 

First measurements of correlated fluctuations from clus- 
tered infrared galaxies at sub-mm wavelengths have been 
detected by (Lagache et al. 2007; Grossan & Smoot 2007; 
Viero et al. 2009; Hall et al. 2010; Dunkley et al. 2010; Hajian 
et al. 2011). Our first results (Amblard et al. 2011) have ex- 
tended these findings at arcminute scales by measuring the 
non-linear 1-halo component for the first time. Modelling 
suggests that at 350 pm 90 per cent of the background in- 
tensity is generated by faint galaxies at 2 > 1 in dark matter 
halos with a minimum mass of log[M m i n /MQ] = 11.5lg 2 , 
agreement with BLAST (Viero et al. 2009). 


3.4 Additional Science Enabled by HerMES 

We expect to detect over 100, 000 sources in our survey. 
The scientific themes explored in sections 3.2 and 3.3 will 
be dramatically extended and improved with the samples 
available now and the full sample once complete. Here we 
mention briefly a very few other science topics that might 
be addressed by us or others using such a large survey. 

The FIR colours of the Herschel sources can help ad- 
dressing the question of how much of the energy produc- 
tion comes from accretion (AGN) and how much from star 
formation. First results on an SDSS sample of AGN (Hatz- 
iminaoglou et al. 2010) find that one third are detected by 
SPIRE, with the long wavelength colours indistinguishable 
from star forming galaxies. Modelling of the full SED re- 
quired the combined contribution of both AGN and star- 
burst components, with the former dominating the emission 
at the MIR wavelengths and the latter contributing mostly 
to the FIR wavelengths. This suggests that SPIRE detects 
the star formation in AGN, with little contamination from 
any dusty torus, offering high hopes for disentangling nu- 
clear and star formation activity. 

The wealth of data in these fields mean we can explore 
the FIR properties of many known samples. Our first results 
on Lyman break galaxies have already shown that we can 
detect U-band dropout sources with stacking (Rigopoulou 
et al. 2010) and FUV drop-out sources individually (Bur- 
garella 2011). We have also shown that galaxies selected on 
the basis of the Spitzer IRAC colours probe a wide range of 
FIR temperatures (Magdis et al. 2010). 

We will compare the FIR measure of star-formation 
with other tracers. In collaboration with the PEP team we 
examined the well-known FIR radio correlation in GOODS- 
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N (Ivison et al. 2010). Exploring giR, i.e. the logarithmic 
ratio of the rest-frame 8-1000 /im flux and the 1.4-GHz flux 
density, there is no evidence that 3 ir changes signicantly for 
the whole sample: qir oc (1 + z)~* , where 7 = — 0.04±0.03 at 
z = 0 — 2, although if the small volume at 2 < 0.5 is removed 
we find 7 = —0.26 ± 0.07. HerMES will create a complete 
data set to understand the global relationship between FIR 
and optical galaxies, the effect of dust attenuation in opti- 
cal/UV populations, and phenomena in individual galaxies. 
First results comparing HerMES and GALEX (Buat et al. 
2010) confirm that total infrared luminosity accounts for 90 
per cent of the total star formation rate, though this reduces 
to 70 per cent when considering the lower star formation rate 
systems (M* < lMgyr -1 ). 

These ancillary data can also be used to investigate the 
detailed properties of the FIR galaxies, e.g. their morphol- 
ogy. One study has explored galaxies with morphological 
classifications at 2 < s < 3 and shows that the mean SFR 
for the spheroidal galaxies is about a factor of three lower 
than for the disk like galaxies (Cava et al. 2010). 

Observations of the rich clusters - the densest known 
regions of the Universe - yield information about their as- 
trophysics and history via the Sunyaev-Zel’dovic (SZ) ef- 
fect (Birkinshaw 1999; Carlstrom et al. 2002), which dom- 
inates the extended several-arcmin-scale emission of clus- 
ters at wavelengths longer than about 500 /im. The SZ effect 
arises from inverse Compton scattering of cosmic microwave 
background photons by hot (1-10 keV) gas in the intraclus- 
ter medium. We intend to combine SPIRE and Planck data 
to measure the SZ effect and the sub- millimetre foregrounds 
between 150 GHz and 1 THz. Based on the different spec- 
tral shapes of the SZ effect and foregrounds, SPIRE data 
will enable us to separate out Galactic dust, cluster and 
background galaxies, the thermal SZ effect and the effects 
of relativistic electrons. 

4 DATA PRODUCTS 

4.1 SPIRE catalogues 

As an illustration of the kind of data products that Her- 
MES will produce we show an approximation to the SPIRE 
250 p m survey areas and depths in Table 7 (together with 
H- ATLAS and GOODS-H). We indicate an estimate of the 
numbers of galaxies on the sky from the Valiante et al. (2009) 
model which is one of the best fits to the current data and 
to a direct determination of the counts from both resolved 
sources (Oliver et al. 2010b) and fluctuation analyses (Glenn 
et al. 2010). Finally we give an estimate of the numbers of 
catalogued sources above those flux density limits estimated 
from our 24 p m driven extractions (at deep levels) and our 
single-band detections at shallow levels. Overall we thus ex- 
pect 100,000 sources detected at > 5a. 

4.2 Ancillary Data 

4-2.1 Required Ancillary data 

To estimate the required ancillary data we have examined 
our first cross-indentified catalogues (Roseboom et al. 2010). 
These are lists with photometry at the positions of known 
24 /im galaxies and thus are not a complete description of 


Levels 

Area 

[deg 2 ] 

5&250 

[mjy] 

Aval. 

[10 3 ] 

AGlenn 

[10 3 ] 

Ac at 
[10 3 ] 

PACS Ul. 

0.012 







Level 1 

0.15 

4 

2.2 

2.0 

± 

0.1 

— 

Levels 2-4 

6.0 

10 

17 

22.4 

± 

0.9 


Level 5 

37 

15 

53 

73.6 

± 

2.3 

52 

Level 6 

52 

26 

20 

28.1 

± 

0.6 

30 

H-ATLAS 

570 

45 

76 

90.6 

± 

2.9 

115 

Level 7 (HeLMS) 

270 

64 

130 
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Table 7. Projected SPIRE survey results for the 250 (im band. 
This table simplifies the survey giving approximate instrumental 
noises in 4 tiers (LI includes GOODS-N). The 5cr confusion noise 
from Nguyen et al. (2010) is 29 mjy, approximately the Level 6 
depth. Numbers of 250 pm sources are estimated from: a count 
model (Valiante et al. 2009, Aval); our P(D) analysis (Glenn 
et al. 2010, Vcienn ) and from our raw number counts in fields 
that we have at these depths, extracted as described in (Smith 
et al. 2011, Meat)- 


Band 

Units 

10% 

50% 

90% 

UV(0.2A) 

[AB] 

22.0 

28.3 

33.7 

R 

[AB] 

18.6 

22.5 

25.2 

I 

[AB] 

18.1 

21.5 

23.6 

K 

[AB] 

17.2 

19.5 

20.8 

3.6 fim 

My] 

380 

90 

30 

24 /im 

My] 

3000 

880 

220 

70 /im 

[mjy] 

42 

13 

4.7 

850 pm 

[mjy] 

6.8 

2.3 

1.1 

21 cm 

My] 

330 

100 

50 


Table 8. Estimates of depth required to detect SPIRE galaxies at 
various other wavelengths. The estimates are based on the mock 
catalogues of Xu et al. (2003) cut to have S 250 > 30mjy. We 
tabulate the depth at which a given percentage of the catalogue 
would be detected. 

the Herschel populations; however they are approximately 
90 per cent complete. 

In Fig. 17 we show the number of sources as a function 
of 250 /im flux and i or K s band magnitude. 

4-2.2 Available Ancillary data 

The survey fields are very well studied and it is outside 
the scope of this paper to provide a complete description 
of all the many ancillary data that are available in these 
fields. A more detailed description of the ancillary data 
will be provided by Vaccari et al. (in prep.). Our inten- 
tion is to homogenise and make publicly available all an- 
cillary/complementary data in our final data release. 

4-2.3 Deliverable data products 

Our intended data products are summarised in Table 9. 
The Herschel source catalogues from SPIRE and PACS data 
(SCAT and PCAT respectively) will consist of the usual 
independent lists where sources are selected from data at 
one wavelength without reference to any other. Associated 
with these catalogues will be validation analyses, including 
completeness, reliability and the information necessary to 
construct selection functions for standard scientific analysis. 
In addition these products will include fluxes estimated for 
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Name 

Description 

Minimum Parameters 

SCAT 

SPIRE Source Catalogues 

Positions, Fluxes, errors, SNRs, etc. 

SMAP 

SPIRE Maps 

Maps of flux, noise and coverage 

PCAT 

PACS Source Catalogues 

Positions, Fluxes, errors, SNRs, etc. 

PMAP 

PACS maps 

Maps of flux, noise and coverage 

SPCAT 

SPIRE/PACS band- merged catalogues 

Positions, Fluxes, errors, SNRs, etc. 

CLUS 

Catalogues &; Maps for Clusters 

As above for maps and catalogues 

XID 

Cross identifications with selected homogenous 
catalogues at other wavelengths. 

Fluxes, errors, SNRs, positions, positional off- 
sets 


Table 9. Deliverable Data Products. 



s 250 (mjy) 


Figure 17. Density of SPIRE sources as a function of 250 /jm flux 
density and optical, /-band (top) and K s (bottom) magnitudes. 
The dashed line indicates the optical or NIR depth required to 
detect 90 per cent of the sample at a given 250 pm flux density, 
while the dot-dashed lines is the depth to detect 10 per cent. 

sources from other catalogues (including sources from public 
Spitzer catalogues). Our first SCAT products are described 
in Smith et al. (2011) and our first PCAT products by Aussel 
et al. (in prep.). 

The SPCAT product will include all Herschel bands. 
Upper limits will be listed for sources detected in some Her- 
schel bands but not others. 

The XID product will include associations with a vari- 
ety of large homogenous catalogues, including, but not nec- 
essarily limited to, public Spitzer catalogues. Our first XID 
products are described by Roseboom et al. (2010). 


Maps from SPIRE and PACS data (SMAP and PMAP 
respectively) will be suitable for extended source analysis, 
fluctuation analysis etc. Our first SMAP products are de- 
scribed by Levenson et al. (2010). 


4-2-4 Other Data Products 

We expect to produce additional data products as an out- 
put of the pursuit of our science goals. These will include 
maps and catalogues of sources from data acquired at other 
facilities (optical, near-IR, radio etc.). It will also include 
value-added products where observed data have been used 
to model other properties of the catalogued objects, such as 
photometric redshift, luminosity or spectral class. It is im- 
possible to define a complete list of such products at this 
stage. We will make these available to the community on a 
best-efforts basis. 


4-2.5 Simulated data 

In order to plan our surveys and simulate our expectations 
we have compiled and homogenised mock catalogues from 
these and other models, which are publicly available via 
hermes . Sussex . ac . uk/. These and other simulations will be 
made available on a best-efforts basis through this site. 


4-2.6 Data Release Schedule 
Early Data Release: EDR 

Our first data release was proposed to be in time for the 
second open call for Herschel proposals (OT2). This was be- 
fore the Science Demonstration Phase (SDP) release rules 
were established and when OT2 was expected to be ear- 
lier. In fact our SDP Early Data Release was made on 
2010 July 1. This meant it was in time for OT1 (due on 

2010 July 22). This data release is described in Smith et al. 
(2011) and, as we proposed, it was restricted to SPIRE high 
signal-to-noise sources in order to be as reliable as possi- 
ble. It included maps from our Abell 2218 observation (# 1) 
and 250 pm catalogues limited at S 250 > 100 mjy for all 
our SDP fields (FLS#40, GOODS-N #14, Lockman-SWIRE 
#28, Lockman-North #19). 

A second Early Data Release EDR2 was made on 

2011 September 19 which included bright source catalogues 
similar to those for EDR but for the DR1 fields (see Table 1.) 
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Data Release 1: DR1 

An extensive Data Release (DR1) of maps and catalogues 
will be made on 2012 March 27. DR1 will include data from 
the SDP observations and all SPIRE observations completed 
by 2010 May 1 (A2219 #7, MS0451.6-0305 #3, ECDFS #15, 
XMM-LSS #36, EGS HerMES #29, Groth Strip #17, 
Bootes #37, ADFS #38, ELAIS N1 HerMES #31). All 
products will be accompanied by documentation in the form 
of papers in refereed journals. 

Data Release 2: DR2 

DR2 will occur at the end of the mission. This will include 
all our deliverable data products and ancillary data in their 
final form. 

4.3 Archival Value and Data Access 

As our observations are in all the most well-studied sur- 
vey fields, the legacy value is enormous. We fully ex- 
pect a rich data-base, leading to abundant science be- 
yond the resources of our team. In addition to any 
ESA data releases (herschel.esac.esa.int/) our data 
will be released through the Herschel Database in Mar- 
seille, HeDaM (hedam.oamp.fr/HerMES). The information 
system design and its implementation are developed under 
the SItools middleware interface provided by the ONES 
(vds . cnes . f r/sitools/). The data (images and catalogues) 
are accessible in various formats (fits files, VOTable, ascii) 
and accessible through Virtual Observatory Tools. Ad- 
vanced searches, cross correlated data and the correspond- 
ing images are also implemented, including visualization fa- 
cilities like ALADIN (http://aladin.u-strasbg.fr/) and 
TOP CAT (http : //www. star .bris . ac .uk/~mbt/topcat/). 


5 DISCUSSION AND CONCLUSION 

We have presented the Herschel Multi-Tiered Extra-galactic 
Survey (HerMES). This survey builds on the legacy of ex- 
isting FIR and sub-mm surveys. It will provide a census of 
star-formation activity over the wavelengths where the ob- 
scured star-formation peaks and over representive volumes 
(and thus environments) of the Universe at different epochs. 
It is being carried out in some of the best studied extra- 
galactic fields on the sky, which is invaluable for the inter- 
pretation of the data both technically, by enabling accurate 
identifications and reducing the impact of confusion noise, 
and scientifically, by allowing exploration of the physical 
processes manifest at different wavelengths. We have pro- 
vided the description and rationale of the survey design. We 
also described the data products we plan to deliver and their 
schedule. 

Our first results from the Science Demonstration Phase 
data have fully demonstrated the promise of the full sur- 
vey. We have quantified the confusion noise at SPIRE wave- 
lengths (Nguyen et al. 2010), 5o250 = 29.0±1.5 injy, finding 
it to be very similar to what was anticipated. This confu- 
sion is challenging to deal with (e.g. Brisbin et al. 2010) but 
we are exploring sophisticated techniques to deal with this 
(e.g. through prior positional information, Roseboom et al. 


2010) and using P(D) analysis have already probed to 4 mjy 
and accounted for 64 per cent of the background at 250 fim 
(Glenn et al. 2010). It seems that previous phenomenolog- 
ical galaxy populations need revision (Oliver et al. 2010b; 
Glenn et al. 2010) and we now anticipate that we will be 
able to catalogue over 100,000 galaxies with > 5<r detections 
at 250 /im. The galaxies appear to be the luminous actively 
star-forming galaxies we expected (e.g. Elbaz et al. 2010) 
with a strongly evolving luminosity function (Vaccari et al. 
2010; Eales et al. 2010b). Also, as expected, SPIRE probes a 
wide range of effective temperatures, including warm galax- 
ies and those cooler galaxies typically seen by sub-mm sur- 
veys (Hwang et al. 2010; Magdis et al. 2010; Chapman et al. 
2010; Roseboom et al. 2011). A clue to the problems that 
the phenomenological models have may lie in the hints of 
the presence of cooler than expected dust in some galax- 
ies (Rowan- Robinson et al. 2010; Schulz et al. 2010). We 
also see evidence for sources being magnified through grav- 
itational lensing by foreground galaxies in the field (Schulz 
et al. 2010; Wang et al. 2011; Conley et al. 2011), and in tar- 
geted clusters. These magnified galaxies provide a window to 
study intrinsically lower luminosity galaxies at higher red- 
shifts. We have identified strong clustering of SPIRE galax- 
ies (e.g. Cooray et al. 2010; Amblard et al. 2011), indicating 
that these luminous systems lie in massive dark matter ha- 
los and implying they are the progenitors of galaxies in rich 
groups and clusters today, i.e. elliptical galaxies. 

HerMES will constitute a lasting legacy to the commu- 
nity, providing an essential complement to multi-wavelength 
surveys in the same fields and providing targets for follow-up 
using many facilities, e.g. ALMA. The results are expected 
to provide an important benchmark for theoretical models 
of galaxy evolution for the foreseeable future. 
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APPENDIX A: DETAILED RATIONAL FOR 
DEFINITION OF EACH SURVEY REGION 

Our deepest tier, Level 1 (#13), covers the GOODS-S region 
which is one of the two deepest Spitzer fields (Dickinson et al. 
2003). 

The other GOODS field, GOODS-N, is covered by one 
of our Level 2 observations (#14), though our observa- 
tions are substantially wider. The boundaries of our other 
Level 2 field, the Extended Chandra Deep Field South field 
(ECDFS, #15), is defined by the deep FIDEL coverage 
(Dickinson & FIDEL team 2007). 

Our Extended Groth Strip (EGS) field at Level 3 
(#17) is also defined to match the FIDEL boundaries. The 
Lockman-East field at Level 3 (#18, #18B) covers Spitzer 
guaranteed time program data (#18) and the Spitzer Legacy 
program of Egami et al. (#18B). Those deep sets (#13, 14, 
15, 17 and 18) were all co-ordinated with the ‘PACS evolu- 
tionary Probe (PEP, Lutz et al. 2011) team. The Lockman- 
North field at Level 3 (#19, 20) covers the deep Spitzer field 
defined e.g. in Owen & Morrison (2008). 

The UDS field at Level 4 (#23) is defined by the Spitzer 
SpUDS observations (Dunlop et al. 2007) and we observe 
this field at Level 3 (#21) with PACS. The Spitzer COSMOS 
field is observed in #22 and # 22B, though our principal 
definition was the PEP observation of this field (discussed 
more in Section 2.8). The VVDS field at Level 4 (#24, 26) 
is not defined by Spitzer observations but by the optical 
spectroscopic survey of Le Fevre et al. (2005). 

The Level 5 and 6 fields CDFS SWIRE, Lock- 
man SWIRE, XMM-LSS SWIRE, ELAIS N1 SWIRE, 
ELAIS N2 SWIRE (#27, 28, 34-36, 39 and 41) are defined by 
the SWIRE fields (Lonsdale et al. 2003) - those fields based 
in turn on the European Large Area ISO Survey, ELAIS, 
(Oliver et al. 2000); the XMM-LSS Survey (Pierre et al. 
2006) and flanking the Chandra Deep Field South (Giac- 
coni et al. 2001) and various Lockman Hole fields (Lockman 
et al. 1986)]. The Bootes NDWFS field at Level 6 (#37) is 
defined by the Spitzer Guaranteed time survey (Jannuzi & 
Dey 1999). The FLS field Level 6 (#40) is defined by the 
Extragalactic part of the Spitzer First Look Survey (Fadda 
et al. 2006) and is commonly referred to now as XFLS. The 
AKARI deep field south (ADFS, #38) is defined with ref- 
erence to the Spitzer (Scott et al. 2010a; Clements et al. 
2011) and BLAST observations (but see Section 2.8). The 
Level 5 observations in #29, 30, 31, 32, 39B lie within or in- 
clude other fields but the bounding regions are new (hence 
labelled ‘HerMES’ or VIDEO) and have been planned with 
the expectation of subsequent follow-up with the SCUBA-2 
Cosmology Legacy Survey (Dunlop et al. 2010), the Spitzer 
SERVS survey (Lacy & SERVS team 2009) and the VISTA- 
VIDEO survey (Bonfield et al. 2010). The fields #29, 32 
and 39B were jointly defined in co-ordination with VISTA- 
VIDEO who fixed the final field location. 


APPENDIX B: MODELLING OF SPIRE 
DITHERING PATTERNS 

SPIRE maps are built by scanning an array of bolometers 
across the sky in a raster with long scan legs each separated 
by a short step, # max (e.g. # m ax = 348" for SPIRE ‘Large 
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Map’ mode). The resulting hit-rate or coverage of detector 
readouts per sky bin is non-uniform (an effect which is ex- 
acerbated by dead or noisy bolometers). This non-uniform 
coverage can be improved by ‘dithering’, i.e. repeating scan 
with offsets. We have modelled this to try and optimise the 
dithering pattern. 

Since we are interested in point sources we can assume 
that the detector readouts will be combined with a point 
source filter (e.g. Smith et al. 2011). The flux estimator for 
a source, / will be given by the 

1 J2i w i d i/ p i 


where di is the readout of detector i, Pi is the point source 
profile for the source at detector i and Wi is a weighting. The 
optimal filter for isolated sources is Wi = Pf/of , where a; is 
the noise of the detector i. The variance in this estimator is 
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We can consider the two scan directions independently 
so we need only model the coverage in one-dimension. The 
sequential scan legs introduce a symmetry on the scale 
flmax, so we project each bolometer position onto the range 
0 < di < flmax in the cross-scan direction. We then construct 
a one-dimensional variance profile V{0) by analogy with 
equation Bl setting Pi the point spread function P(d — di). 
For simplicity we set Wi = 0 for dead or noise bolome- 
ters and <Ji = 1 otherwise and used Gaussian beams with 
FWHM=18. 15/25. 15/36. 3” for the three bands. 

We then defined a metric, M, to optimise dither pat- 
terns, on the understanding that we want to reduce the vari- 
ation in variance. When considering the dither pattern for 
one band in isolation we simply used the fractional variance 
of the variance 



where the sum is over the profile. As the SPIRE bolome- 
ters scan the sky simultaneously in all bands any dithering 
scheme would apply to all bands. However, considering three 
bands simultaneously there is no obvious metric (unless we 
considered sources of a particular colour) we did define an 
arbitrary metric M 2 = Mp SW + Mp MW + Mp LW but have 
restricted this discussion to the single bands independently. 

The aim is to choose a dither pattern that minimises 
M. If N independent scan maps with N—l dither positions 
are available then the dither pattern is defined by N — 1 
offsets A d. We adopted four approaches: (a) optimisation 
by brute-force search through N—l dimensional space (only 
attempted up to N = 4) (b) sequential optimisation where 
we chose the best Adi for each additional dither, i, given 
the A 9 found for the previous i — 1 dithers (c) equal spacing 
Adi = A #2 = • • . = d max /N (d) random spacing with Adi 
uniformly selected from 0 < Adi < d ma , x . 

For low values of N < 4 where both were calculated 
we found that the brute-force optimisation (a) agreed rea- 
sonably well with the sequential optimisation (b). We found 
that the equal spacing (c) performed similarly to the se- 
quential optimisation at low N and typically better at high 
N > 10) except at specific N (e.g. N = 15 for PSW and 


#max = 348”) when the projected bolometer spacing were in 
phase. Random offsets (d) were invariably worst. The raw 
variation with no dithers (N = 1) was 12, 15, 10 per cent for 
PSW, PMW and PLW respectively, this declined rapidly to 
about 3 per cent by N = 3 and was < 1 per cent for N > 16. 

A penalty for dithering with these large steps is that 
the ramp down in coverage at the edges of the map is more 
gradual, i.e. less area at the full coverage with more area at 
low coverage. When designing offsets in both scan directions 
we chose pairs of offsets tracing a square to reduce the im- 
pact of this ramp-down and this strategy is included in the 
SPIRE Observers’ Manual. 



